Mechanism of membrane tube formation induced by adhesive nanocomponents 
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We report numerical simulations of membrane tubulation driven by large colloidal particles. Using 
Monte Carlo simulations we study how the process depends on particle size, concentration and 
binding strength, and present accurate free energy calculations to sort out how tube formation 
compares with the competing budding process. We find that tube formation is a result of the 
collective behavior of the particles adhering on the surface, and it occurs for binding strengths that 
are smaller than those required for budding. We also find that long linear aggregates of particles 
forming on the membrane surface act as nucleation seeds for tubulation by lowering the free energy 
barrier associated to the process. 



A key factor in cell trafficking and intercellular com- 
munication is the internalization of complex macro- 
molecules. As large and charged biological cargo cannot 
directly cross the lipid bilayer that envelops the differ- 
ent compartments within eukaryotic cells, this process 
is usually accompanied by the formation of vesicular- 
and tubular-shaped membrane protrusions. The mech- 
anism by which they develop can be extremely diverse 
[T]-[6]. It often involves active processes requiring acces- 
sory factors, such as clathrin or caveolin protein coats, 
or motor-proteins and external forces. It can also de- 
velop as a result of the self-assembly of anchoring pro- 
teins, such as BAR domain proteins [3(8], that impose a 
local curvature on the lipid bilayer. The physical mech- 
anism driving protein aggregation in this case is fairly 
well understood within the framework of effective bend- 
ing mediated Casimir forces [9 j. The size and shape of the 
resulting deformation is determined by how the packing 
properties of the proteins couple to the elastic response 
of the membrane. 

Several endocytic pathways, however, are found to be 
triggered by the cargo itself [2j HJ HOI E] • In some cases, 
such as HIV-1 [12] , the virus itself is formed on the mem- 
brane as its proteins self-assemble inducing their own 
vesicular bud. The internalization is thus a consequence 
of cooperativity of many protein molecules. In this paper 
we are interested in passive internalization of preassem- 
bled viruses, virus-like particles and other colloidal parti- 
cles. The main difference from the cases discussed above 
is that the interaction of a single colloidal particle (typi- 
cally one order of magnitude larger than a protein) with 
the lipid bilayer can induce its own invagination by wrap- 
ping its surface with the membrane. For instance, it has 
been shown that budding of preassembled alphaviruses 
and type-D retroviruses [U [13] , as well as charged col- 
loids [14] , can take place without the presence of external 
factors. 

Although one might expect budding to be the main 
mechanism for internalization of large particles, long 
tubular protrusions typically of one-particle diameter are 
often observed in viral or nanoparticle internalization. 
Simian virus 40 (SV40), upon its binding to membrane 
receptors, is found to induce deep invagination and tubu- 



lation of both the plasma membrane and giant unilamel- 
lar vesicles (GUVs) [15] . Its entry occurs via small, tight- 
fitting indentations and the resulting invaginations have 
the same size as the virus-particle diameter. Positively 
charged nanoparticles were also shown to spontaneously 
induce tubulation in supported [16] and unsupported 
[T4] giant unilamellar vesicles, suggesting the existence of 
a general physical mechanism of internalization, which is 
not exclusive for viruses and does not require assistance 
of membrane proteins. 

Understanding this phenomenon is of great importance 
for developing ant i- viral strategies, but also because vi- 
ral and virus-like particles, as well as artificial nanopar- 
ticles, are promising tools in gene-therapy and molecular 
medicine, for which control over their cellular uptake is 
essential. Despite the large body of work [T7H24] on the 
particle budding problem, most studies have focused on 
the interaction of a single particle with the membrane, 
and have completely missed tube formation, a crucial 
component of the phenomenological behavior associated 
to particle internalization, that can only arise as a re- 
sult of nontrivial cooperative behavior among many par- 
ticles. Here we use computer simulations to investigate 
the physical mechanisms behind the occurrence of this 
process, and show how it depends on particle size, con- 
centration and binding strength. While the phenomenon 
has been observed in several experiments, to the best 
of our knowledge, this paper presents the first theoreti- 
cal study that addresses nanoparticle-driven tubulation, 
and rationalizes its interplay with the particle budding 
process. 

Our system setup consists of N p particles, modeling 
colloidal viruses, virus-like particles or inorganic colloids, 
placed inside a vesicle of undeformed average radius R. 
Given the large size difference between the thickness of 
the vesicle and the nanoparticles considered in this study, 
we model the vesicle using a simple one particle-thick 
solvent-free model consisting of N spherical beads of di- 
ameter (7 connected by flexible links to form a triangu- 
lated network [25H27] whose connectivity is dynamically 
rearranged to simulate the fluidity of the membrane. The 
membrane bending energy acts on neighboring triangles, 
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and has the form 
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where ^ is the bending rigidity, and the and nj are 
the normals of two triangles i and j sharing a common 
edge. The cost associated with area changes is included 
via the energy term = jA, where 7 is the tension of 
the surface and A is its total area. The particles are rep- 
resented as spheres of diameter cr np = Zcr, where Z > 1 
is a parameter used to control their size. Excluded vol- 
ume between any two spheres in the system (particles 
and surface beads) is enforced with a hard-sphere poten- 
tial. Finally, the colloid-to-membrane adhesion energy is 
modeled via an additional power-law interaction between 
the particles and the surface beads defined as 
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and cut-off at r cut = 1.5a M , with gm = (cr -\- cr np )/2. Do 
is thus the membrane-particle binding constant. This 
potential is quite generic and is typically employed to 
describe short-range attractions, such as ligand-receptor 
or van der Waals interactions. The system is equilibrated 
using the Monte Carlo simulations in the NVT ensemble 
and most of our data are obtained at ^5 = 5&bT, 
7 = lk B T/a 2 and Z = 2, 3, 4, 6 or 8. 

To qualitatively understand the interaction between a 
single particle and a membrane, consider a particle of 
radius R p having a spherical cap of height h and area 
Scap = 2irR p h in contact with a membrane. The elastic 

costs associated with this configuration are ^S cap and 

p 

77r/i 2 , where the first term accounts for the bending and 
the second for the surface energy of the membrane. The 
free energy gain due to the binding energy between parti- 
cle and membrane scales as — DoS CSiP . A balance of these 
terms leads to an equilibrium particle-membrane inden- 
tation, /i, and to a particle coverage xt = S cap / (4ttR p ) = 
D fcb/(2i? p ) ^ rpj^-g SU gg es tg that ^he particle will bud off 

the membrane as soon as D > k h /(2R 2 p ). Although more 
sophisticated calculations have been put forward to un- 
derstand the nature of this transition [I8j [28-30 , the 
main result is that for a given binding constant Do, bud- 
ding is easier for large particles. This scaling argument 
gives a simple explanation of why this process is likely for 
colloidal particles and preassembled viruses, but not for 
single proteins and small macromolecules, and provides 
us with an intuitive framework from which to begin our 
study. 

We begin our analysis by computing a diagram that in- 
dicates, for a given value of Do an d R pi the phenomeno- 
logical behavior of the particle-membrane coupled sys- 
tem at a constant particle concentration. The results 
are shown in Fig. [I] For small values of Do, the overall 
shape of the membrane is unchanged while the particles, 
barely adhering to it, freely diffuse over its surface as a 



low-density two-dimensional gas (G). Increasing Do, we 
find that the nanoparticles organize into linear aggregates 
(L). This phase develops due to effective interactions be- 
tween the particles driven by the membrane's need to 
minimize its elastic energy while maximizing its bind- 
ing surface to the particles (see |3lJ [32] for a detailed 
analysis of this phase and its experimental observation). 
Upon further increase in Do, spontaneous formation of 
tubular protrusions (T) takes place. This region of the 
diagram is characterized by nanoparticles tightly and lin- 
early packed into tubular structures extruding out of the 
membrane core. The radius of the tubes equals the diam- 
eter of the particles. This behavior is in agreement with 
the SV40-induced membrane invaginations [15], where 
one-particle-wide tubes were also observed, but tubu- 
lation failed to occur if the adhering viruses were un- 
able to form a sufficient number of interactions with the 
membrane binding sites. Further increase in Do causes 
nanoparticles to promptly adhere to the membrane and 
become completely enveloped into a bud before any sig- 
nificant particle diffusion can occur. The T-B transition 
is not abrupt, and a mixture of both "corrugated" tubes 
and single-particle buds is found in the borderline area 
between the two phases. Although in our model buds 
cannot physically detach from the membrane, they are 
easily identifiable by their complete surface-coverage and 
the characteristic sharp membrane neck shape. A single 
particle bud is shown in the inset of Fig. [T] 

The most important message arising from our analysis 
is that tubulation develops as a result of the interaction 
of many particles and should be expected for interme- 
diate binding constants. Such a behavior occurs for all 
particle sizes considered in this study, indicating what 
sets the tube size is the particle diameter and not the 
natural length-sc ale assoc iated to membrane tube for- 
mation, Ro = \J ^5/(27), obtained by pulling experi- 
ments [33, 34 . Moreover, preassembly of nanoparticles 
into linear aggregates seems to greatly facilitate the for- 
mation of long tubes. 

To obtain more physical insight into the mechanism by 
which tubular protrusions form, we considered a series 
free energy calculations. First, we measure the effective 
interaction between two colloids adhering to the mem- 
brane in the T-region of the phase diagram. Using the 
umbrella sampling method [35] , we compute the proba- 
bility that the two particles are at any given separation 
d from each other and estimate the free energy difference 
AF = F(d) - F(oo). Fig. |2(a)| shows AF as a function 
of d, while the inset monitors how the orientation of the 
dimer with respect to the membrane surface, (/?, depends 
on the same variable. This result is quite revealing; the 
elastic cost required to bring together two large mem- 
brane deformations, responsible for the weak mid-range 
repulsion, is replaced by a large energy gain when the 
particles are in contact. The corresponding configura- 
tion is characterized by two particles contained within 
a membrane tube oriented perpendicularly to the mem- 
brane surface. As we have not imposed any constraint on 
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the values of cp, this is clear evidence, at least at the two- 
particle level, that in this region of the phase diagram, 
tube formation is more favorable than budding. 

Using the same procedure, we can also measure the 
free energy as a function of separation between a two- 
particle-tube and a third isolated particle. Our data, 
shown in Fig. |2(b)| tells us that the lowest free energy is 
again achieved when the three particles are in contact in 
a tubular formation. This very important result indicates 
that tubes and free particles bound to the membrane at- 
tract each other, and once a tube is formed, its growth 
by particle addition drives the system towards a lower 
free energy. In both cases, the extent of the repulsion 
and attraction is dependent on the specific region of the 
phase diagram they are computed at. The characteris- 
tic energy barrier at mid-range distance becomes more 
significant as Do increases and the system crosses over 
to the budding regime, implying that for large Do parti- 
cle aggregation becomes rare, making budding the most 
likely barrier-crossing mechanism. This is a kinetically 
dominated regime: in fact, once the budding threshold 
has been overcome, particles would leave the membrane 
before having the time to aggregate. 

Interestingly, in most of our simulations in the T phase 
we observe that tube formation if often preceeded, in par- 
ticular at higher particle densities, by the formation of 
long linear aggregates that eventually extrude from the 
membrane via a tilting mechanism illustrated in Fig. [3] 
This two-step process becomes more significant as we 
move closer to the L-T boundary, suggesting that these 
aggregates function as nucleation seeds promoting the 
transition. To support this idea we perform two sets of 
simulations: in the first set we start from a connected 
four-particle- long linear aggregate, and measure its sur- 
face coverage xt as a function of Do until a tube is 
formed, in the second set we start from an already tubu- 
lated structure and we decrease Do until the tubule re- 
tracts. As shown in Fig. [3J tubulation is accompanied 



by a sudden jump in the particle coverage Xt (and con- 
sequentially in the binding energy), indicating the pres- 
ence of a free energy barrier between the two states that 
needs to be crossed for the linear aggregates to protrude 
out of the membrane. This result is consistent with 
previous force-extension calculations and experiments on 
GUVs [33], that also indicated tube formation to be a 
first order transition. Finally, we measured the onset 
value Dq at which a preformed linear aggregate forms a 
tube as a function of its size, at a fixed particle radius. A 
weak but clearly inverse dependency is found, shown in 
Fig. |1J and supports the idea that the free energy cost for 
tubulation from the L phase does indeed decrease mono- 
tonically with the size of the aggregates which therefore 
act as nucleation seeds for the transition. It should be 
stressed that because the probability of forming linear ag- 
gregates increases with particle concentration, it is logical 
to expect tubulation to be more likely to occur in denser 
systems. This is indeed what we find in our study. 

We have shown that for a broad range of binding en- 
ergies, tube formation and not membrane budding is the 
main mechanism leading to internalization of sufficiently 
large particles. Nowhere in our simulations have we ob- 
served formation of membrane tubes of radius larger than 
one particle diameter; however, these may develop as a 
result of direct particle-particle interactions or nontrivial 
long-range electrostatic effects [14] not included in our 
study. It should be emphasized that our results should 
hold as long as the particle size is sufficiently large so 
that the molecular details of the membrane can be ig- 
nored. Although the elastic constants of our model were 
selected in a range relevant to biological processes and we 
only considered two vesicle radii, we do not expect the 
process to be extremely sensitive to these parameters. It 
would be useful to carry out a surface minimization anal- 
ysis of this problem as it was performed for force driven 
tubulation [34], and we really hope that our work will 
stimulate further experimental work in this direction. 
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FIG. 1: Left panel: Do-R p phase diagram of the membrane aggregates and protrusions induced by colloidal particles. Right 
panel: Snapshots of the linear (L) and tubular (T) phases. The inset shows a typical single-particle bud conformation (B) that 
occurs at large Do- The bottom region of the phase diagram is the gaseous phase (G). The radius of the membrane is R — 30a 
and the particle surface fraction is kept constant at 0.15. 




FIG. 2: Tube formation and growth, (a) Free-energy as a function of separation of two membrane-bound particles. The inset 
shows the orientation of the dimer with respect to the membrane surface and indicates the distance at which the tubulation 
occurs. Here, ip is defined as the angle formed by the vector connecting the centers of the two particles and that connecting 
the center of the membrane to the particles' midpoint, (b) Free-energy as a function of separation of a two-particle tube and 
a single membrane-bound particle. In both cases Rp = 4, R = 15a and Do = 2.6/cbT. 
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D [k B 7] 

FIG. 3: Hysteresis associated with the tubulation of a linear aggregate, in terms of the surface coverage \t and Do, for the 
extrusion of a four-particle-long aggregate. \t is computed as the ratio between the number of membrane beads in contact 
with the particles and the same number when the surface completely envelops the particles. The red crosses show the results of 
simulations that start form a linear aggregate, while the black circles show simulations that start from a tube. Here Rp = 4a 
and R — 15a 
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FIG. 4: Onset binding constant for tubulation Dq as a function of the length of the preformed linear nucleation cluster; Rp = 4cr 
and R — 15a. 



